Multidimensional cut-off technique, 
odd-dimensional Epstein zeta functions and 
Casimir energy of massless scalar fields 



*Ariel Edery 
Bishop's University, Physics Department 
Lennoxville, Quebec 
J1M1Z7 



Abstract 

Quantum fluctuations of massless scalar fields represented by quantum 
fluctuations of the quasiparticle vacuum in a zero-temperature dilute 
Bose-Einstein condensate may well provide the first experimental arena 
for measuring the Casimir force of a field other than the electromag- 
netic field. This would constitute a real Casimir force measurement 
- due to quantum fluctuations - in contrast to thermal fluctuation ef- 
fects. We develop a multidimensional cut-off technique for calculating 
the Casimir energy of massless scalar fields in d-dimensional rectangu- 
lar spaces with q large dimensions and d—q dimensions of length L and 
generalize the technique to arbitrary lengths. We explicitly evaluate 
the multidimensional remainder and express it in a form that converges 
exponentially fast. Together with the compact analytical formulas we 
derive, the numerical results are exact and easy to obtain. Most impor- 
tantly we show that the division between analytical and remainder is 
not arbitrary but has a natural physical interpretation. The analytical 
part can be viewed as the sum of individual parallel plate energies and 
the remainder as an interaction energy. In a separate procedure, via 
results from number theory, we express some odd-dimensional homo- 
geneous Epstein zeta functions as products of one-dimensional sums 
plus a tiny remainder and calculate from them the Casimir energy via 
zeta function regularization. 
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1 Introduction 



The Casimir force remained for a long time one of the more esoteric forces in 
Physics attracting at best some theoretical interest. All of this has changed 
in the last eight years or so. After nearly 50 years since its prediction in 1948 
by Casimir pP, the force has now been successfully measured by a modern 
series of experiments starting with Lamoreaux's 1997 landmark experiment 
[3J with a torsion pendulum which reduced errors dramatically compared 
to the early 1958 experiment by Spaarnay 2 . The force was subsequently 
measured more precisely in 1998 using an atomic force microscope [I] and 
the measurements agreed with theoretical predictions to within 1% after 
finite conductivity, roughness and temperature corrections were taken into 
account. Thus the modern era of precise Casimir measurements was born 
and a non-exhaustive list of other experimental studies since then can be 
found inJHlliBiaillllBinjIIlQS]. Interest in the Casimir force has also 
been fueled by theories with large extra dimensions which predict among 
other things a deviation from Newtonian gravitation at the sub-millimeter 
scale |14j . To date no deviation has been found. Recently, a Casimir force 
experiment ^H] has placed new constraints on the parameters of such pro- 
posed theories. An up-to-date list of gravitational experiments can also be 
found in 15 . As with many fundamental Physics discoveries, at first the 
Casimir force seemed to have no apparent engineering application (since it 
is significant only on micron or submicron scales). However, our ever in- 
creasing ability to build structures on smaller scales has made the Casimir 
force something various industries need to take into account. For example, 
in 2001, scientists at Lucent Technologies showed that the Casimir force 
could be used to control the mechanical motion of a microelectromechanical 
system (MEMS) device ^7j( see a ^ so the recent paper ^1 and references 
therein). MEMS are micron-sized devices in which tiny sensors and actua- 
tors are carved into a silicon substrate and are currently in use as car air-bag 
sensors. For more details on the Casimir effect the reader is referred to the 
following books [13 H3 and reviews [301 I3T | 132 1 I33| . 

All the measurements of the Casimir force to date have been limited to the 
case of the electromagnetic field. However, experiments may soon (or may 
already have done so indirectly) measure the Casimir force for a massless 
scalar field. Quantum fluctuations of the quasiparticle vacuum in a zero- 
temperature dilute Bose-Einstein condensate (BEC) should give rise to a 
measurable Casimir force as explained in recent papers |18l I19j . The au- 
thors in ^HJ QUI state that indirect effects from these quantum fluctuations 
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may have already been observed |2L)l 1211 1221 127j . Note that this is a real 
Casimir effect due to quantum fluctuations in contrast to thermal fluctua- 
tions (often called pseudo- Casimir). The fact that the field propagates at 
the speed of sound in the BEC medium in contrast to the speed of light 
in Minkowski spacetime does not change anything fundamental in relation 
to the Casimir energy. If the speed of propagation is constant in a given 
medium, the Casimir energy in units of this speed will be the same value 
regardless of whether the medium is spacetime or a BEC. Moreover, a gen- 
erally covariant action analogous to what we see in General Relativity exists 
for scalar fields propagating in a particular fluid. The Lagrangian is similar 
to that of a massless Klein-Gordon field with the Minkowski metric 77^ v of 
spacetime replaced by an effective or acoustic metric g^y [34] . Quoting di- 
rectly from [HE], "at low momenta linearized excitations of the phase of the 
condensate wavefunction obey a (3+l)-dimensional d'Alembertian equation 
coupling to a (3+l)-dimensional Lorentzian-signature 'effective metric' that 
is generic, and depends algebraically on the background field.". In ^F3\ the 
authors make the important observation that though the dispersion relation 
for quantum fluctuations in a BEC is nonlinear, the Casimir energy picks 
out mostly the long wavelength linear behaviour. This is why the Casimir 
force Fb ec calculated by the same authors |18| I19j for infinitely thin and in- 
finitely repulsive plates immersed in a zero-temperature three-dimensional 
dilute condensate turns out to leading order to be the same as that of a 
massless scalar field moving with the speed of sound v. 

In this paper we are interested in the Casimir effect of massless scalar fields 
traveling with speed v in rectangular cavities of d spatial dimensions where 
q dimensions are large and d — q dimensions are of equal length L. The 
case of arbitrary lengths is also considered in appendix B. We develop a 
multidimensional cut-off technique to solve this problem. Why use a cut- 
off technique? Clearly, it is less efficient than the zeta function technique 
that yields quickly, via analytic continuation, finite results for rectangular 
cavities in terms of Epstein zeta functions. There are a few reasons for the 
importance of the exponential cut-off technique. First, it remains the most 
physically intuitive method. For this reason, recent texts in String Theory 
or Quantum Field Theory (QFT) as well as courses in QFT introduce the 
standard Casimir energy calculation of a string or parallel plates using an 
exponential cut-off. For example, in the text String Theory, Vol. I |23j . 
the Casimir energy for the Bosonic string is handled with an exponential 
cut-off. The result 5Z^=i n ~ > T5" ^ s obtained by replacing n by ne~ Xn and 
extracting the finite result from the series — A + 0{\ 2 ). This cut-off 



3 



method was used instead of the zeta function technique which yields quickly 
£(— 1) = j^. In his recent book, Quantum Field Theory in a Nutshell |24j . 
Zee brings in some humour in explaining a Physicist's perspective on the 
same sum. I quote from p. 66, "Aagh! What do we do with Y^Li n ? None 
of the ancient Greeks from Zeno on could tell us. What they should tell us 
is that we are doing Physics... Physical plates cannot keep arbitrarily high 
frequencies from leaking out.". He then introduces the exponential cut-off 
to damp the ultraviolet frequencies. In the classic QFT text by Itzykson and 
Zuber the electromagnetic parallel plate problem in three dimensions is 
solved via a cut-off function and the Euler-Maclaurin formula and the same 
technique can be seen applied in recent graduate courses (e.g. see "Rela- 
tivistic Quantum Field Theory I, Spring 2003" 26 ). Physicists are therefore 
likely to be familiar with the cut-off technique. Secondly, a multidimensional 
cut-off calculation with an exact determination of the multidimensional re- 
mainder term does not seem to have been systematically carried out for rect- 
angular cavities in arbitrary d dimensions. Papers on Casimir energies in 
arbitrary d dimensions in rectangular cavities have made use of dimensional 
and zeta function regularization |42l 143 \ Hl| . Explicit formulas using the 
exponential cut-off technique in rectangular cavities include parallel plates 
in higher dimensions [3E], rectangular cavities in two and three dimensions 
[S3 EE1 HUE EDI , and explicit formulas via Poisson's formula up to d = 2 ap- 
pear in In |36[ I37[ l4*T] the connection between cut-off and zeta function 
technique is also elaborated and explained. A detailed numerical analysis 
for the electromagnetic case in three-dimensional rectangular cavities can be 
found in Last but not least, by applying the cut-off technique to rect- 
angular cavities we are led in a natural fashion to excellent finite analytical 
formulas plus a remainder. We show that the division between analytical 
and remainder is not some ad-hoc division. The analytical part has a clear 
physical interpretation as sums of parallel plates out of which the rectan- 
gular cavity is constructed. Moreover, the numerical results are excellent 
because the analytical part is trivial to evaluate and the multidimensional 
remainder is derived in a form that converges quickly (exponentially fast). 
As already mentioned, the zeta function technique applied to rectangular 
spaces has the great advantage of leading quickly to finite results expressed 
in terms of Epstein zeta functions. However, one then needs to go a few 
steps further if one wants to express these in a convenient analytical form 
and this is usually a separate procedure. In contrast, analytical results are 
often a natural spin-off of the cut-off technique. 

One section of this paper is devoted to developing a technique that derives 
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highly accurate analytical formulas for a few odd-dimensional homogeneous 
Epstein zeta functions. It turns out that in even dimensions less than or 
equal to 8 one can obtain compact analytical expressions for the homoge- 
neous Epstein zeta function purely in terms of products of one-dimensional 
sums. There is no remainder for these cases. This can be accomplished 
via number theoretic formulas for the representation of integers as a sum of 
squares in even dimensions. For even dimensions above 8, the number the- 
oretic formulas get more complicated and in odd dimensions above 7 they 
are not presently known. For 3, 5 and 7 dimensions the number theoretic 
formulas have only recently been found |50| but they are much more com- 
plicated than in even dimensions. We therefore develop a procedure that 
uses the exact even-dimensional results from number theory and then ap- 
ply the Euler-Maclaurin formula to obtain the odd dimensions. This yields 
the homogeneous Epstein zeta function in 3, 5 and 7 dimensions as a finite 
number of products of one-dimensional sums plus a small remainder term. 
This remainder is even smaller than the remainder obtained via our multi- 
dimensional cut-off technique. For the most important case of 3 dimensions, 
we obtain both a highly compact and extremely accurate analytical expres- 
sion that contains only four terms and where the remainder is a negligible 
0.04% of the Casimir energy. Our specific procedure leads to low remainders 
but is limited to a few homogeneous Epstein zeta functions, albeit one that 
includes the three-dimensional case. A different more general procedure ap- 
plicable to any multidimensional inhomogeneous Epstein-type zeta function 
can be found in |46j . 

2 Multidimensional cut-off technique including re- 
mainder 

In this section we develop a multidimensional cut-off technique to obtain 
formulas for the Casimir energy of a massless scalar field (j>{x) moving with 
a wave velocity v in a d-dimensional rectangular cavity with d — q sides 
of equal length L and q sides of much larger length L m >> L where m 
runs from 1 to q. One can generalize our method to arbitrary lengths and 
this is done in appendix B. Here and throughout the paper we consider the 
more special case as it makes the method, the formulas and the physical 
interpretation more transparent. This section and appendix A (where the 
remainder is evaluated) go together. 
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We consider periodic, Neumann and Dirichlet boundary conditions. The 
fields are assumed to propagate in a homogeneous medium with a constant 
speed v and with a wavelength long enough that the dispersion relation is 
linear i.e. to = vk where k is the wavenumber. In other words, we assume 
the scalar field (j>(x) to obey the standard linear wave equation: 

^>-» 2 WM = o. (i) 

The boundary conditions are either periodic, 4>(x % = 0) = 4>(x % = L), Neu- 
mann, d l <p(x) = at x l = and x l = L or Dirichlet (f>(x l = 0) = 4>(x l = L) = 0. 
Here i runs from 1 to d — q inclusively. After the standard fourier decom- 
position one obtains the following quantized frequencies ui for periodic (p), 
Neumann (N) and Dirichlet (D) conditions : 

Wp = 2 ™(£! + ... + ^ + %!±i + . 
^,D = ™(£l + --- + % i + ! %f ±1 + 

where the rij's run from — oo to oo for periodic boundary conditions, to 
oo for Neumann and 1 to oo for Dirichlet. From quantum field theory we 
know that after quantization the vacuum energy is given by the sum over 
all modes of \ u> (we work in units where ft = 1). The vacuum energies E vac 
for the three boundary conditions labeled (p,N,D) are therefore: 

OO 2 2 

E P vac = ^ E K + --- + -i, + ^fr L + - 

Tli = — OO 

i=l,...,d 
oo 

F vac_ 7TV ST (n 2 -I kn 2 I n d- q +i L 

N ' D ~ 2L ^ ^ 1 d_<? 

rti=0,l 
i=l,...,d 

The above sums are ultraviolet divergent and require regularization. There 
are many different regularization schemes such as exponential cut-off, zeta 
function and dimensional regularization. In this paper the goal is to develop 
a multidimensional cut-off technique via the Euler-Maclaurin formula. Via 
this technique, we obtain formulas for the Casimir energy as a finite sum 
over analytical terms plus a remainder. We fully evaluate the remainder 
term and express it as sums over Bessel functions. We later show that the 
analytical part has an intuitive physical picture: it is the energy needed to 
construct the rectangular cavity out of adding successive parallel plates. We 
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begin by calculating the regularized vacuum energy for periodic boundary 
conditions. After regularization, we then extract the finite Casimir energy 
E p which is the difference between the regularized energy with boundaries 
(discrete modes) minus the regularized energy without boundaries (contin- 
uous modes). We later compare E p to the Epstein zeta function obtained 
via zeta function regularization. We can express the Neumann and Dirich- 
let energies, En and Ed, in terms of sums over E p so only the periodic 
case needs to be evaluated fully. The regularized vacuum energy E reg for 
periodic boundary conditions using an exponential cut-off is: 

OO 11 

rii=—oo 
i=l„„,d 



»?+-+»3-,+^-+-+-fe- (4) 



^n^£ /_ 



-oo 
--l,...,d-q 



oo 



where we replaced the sums over the q large dimensions by integration. The 
parameter A is a free parameter which we later set to 0. The goal is to 
evaluate the expression in (@J) that includes d — q sums and q integrals. Our 
procedure will be to express as an expansion over a function A and 
then use the Euler-Maclaurin formula to evaluate this function. Define the 
following short-hand form for a j — q dimensional sum over q integrals: 



-oo 



oo 



=i,— J-g 



where j runs from q to d — 1 (the case j = q corresponds to no sums, only q 
integrals) . The reader may wonder why we chose a definition with j — q sums 
instead of just simply j. The reason is that the total number of sums plus 
integrals is then j and this simplifies things later on. We define a function 
A by adding one more sum to the above definition: 

rq 



'OO 



£ / e- A dxi . . . dXq 



(6) 
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where the last sum over n excludes zero. With these definitions, we make 
the following useful expansion of 



rii=— oo 
i=l,...,d—q 



+ e7' + e'e/' + e' e 2 /' + - + e e^ 1 



3=1 3=1 

(7) 

Substituting Q into © yields the regularized Casimir energy: 

E^(q,\) = --^flL l U | 9 + ^a A A,(g,A)J . (8) 

i=l y j=g J 

In the above expression, we need to separate the divergent part due to the 
continuum from the finite part related to the Casimir energy as A— >0. The 
term d\ f q contains no sums, only multiple integrals. It is immediately clear 
that this term contributes purely a continuum divergent part as A — > and 
hence makes no contribution to the finite Casimir energy. We now need to 
find an expression for Aj(q, A) given by (jHJ) and extract the finite part related 
to it. To this end we apply the Euler-Maclaurin formula that converts sums 
to integrals. The Euler-Maclaurin formula is given by 51 : 



/ f{x)dx--m-Y,TT\ T i? 2 P / (2p - 1) (0) + J R s 

n=l Jo 2 p=l [2P)1 



(9) 



where /(^"^(O) are odd derivatives evaluated at zero and s is a positive 
integer. The form above for the Euler-Maclaurin formula assumes that the 
function fin) and its derivatives are zero at infinity. R s is the remainder 
term given by |51j 

„1 oo 

where B2 S {x) are Bernoulli functions and f 2s (x + v) are even derivatives of 
/ with respect to x. 
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In applying the Euler-Maclaurin formula to determine Aj(q, A), the function 
/ in question is the exponential function appearing in ©. Regardless of 
the value of p, this exponential function has the property that / 2p_1 (0) is 
zero for all sums in © except the last one over n. A proof of this is given 
in the appendix of |53| . If / 2p_1 (0) is zero for all p it follows that the sum 
from p = 1 to s in @ is zero independent of s. This implies that R s given 
by (jl()j) has the same value for any given s for the case of our exponential 
function. This is proven explicitly in the appendix of |54j . For calculations 
we can simply choose s equal to 1. Since / 2p_1 (0) is zero for all sums except 
the last one, the Euler-Maclaurin formula for those sums reduces to 

00 roc 1 If' 1 00 d 2 

£/(")=/ f(x)dx-~f(0)-~ £ 2 (x)^— f(x + v)dx (11) 

71=1 J ° J ° V=0 

where B2(x) = x 2 — x + 1/6. The function / in © has the property /(rej) = 
/(— rij). The sum over a given nj can therefore be written as 

00 00 

f(TH) = 2 J2 + /(0) 



COO I 

2 I IJ(x)dx--f(0) 



2 1 B2 ^ E h f[x+v)dx ) +/(0) (12) 



/oo 
f(x) dx-R 
-oo 

where R is a remainder given by 

i?= / B 2 (x) -r^f(x + v)dx. (13) 
Jo u=0 dx 

From H12|) we see that each sum in @, except the last one, can be replaced 
by an integral minus R. We therefore have the operator prescription ^ — ► 
J —R. Applying the operator j — q times and then inserting the result in 
© yields 




ir=(/-*) =/ +EM>i\7)r ' *■ <") 



j-q-m 
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and 



21+1 E / e ~ A V n2+ ^ + - +:E ? dan ■ ■ ■ dxj + Rj(q, A) 

n=l ^ 



(15) 



where Rj(q,X) is a remainder given by 

%A) = EE(-l) B 2 J m ? ^- (16) 

m=in=l ^ ' ^ 

Substituting 72 given by lfT5)l into (jTH|l yields 

J— g 00 . f .i m 00 „ 2 

%»=ESHr(';')^'/ / nE/*4 

m=ln=l JU JU i=lvi=0 1 (17) 

e- A > /» a +(* 1 +^) a +-+(*™+^) a +^+--^- m toy... dx m dy x ... d Vj . m 

where the integrations from to 1 and to oo are over the x's and y's 
respectively. The function A given by 1)15(1 contains two terms. The first 
term leads to the analytical part and the second term Rj(q,X) yields the 
remainder. In the limit A = 0, Rj(q,X) is zero but not its derivative with 
respect to A. It is the derivative with respect to A that enters into the Casimir 
energy There is therefore a non-zero contribution to the Casimir energy 
coming from the remainder term and we fully evaluate it later on. For now, 
let us evaluate the analytical term in (|15[). It can be reduced to an infinite 
sum over the modified Bessel function Ko(Xn) which has a useful series 
expansion. We first note that the integral in l)15|) can be expressed in terms 
of the modified Bessel function Kj-i (A n) |52| : 



e - x y n2+xl+ - +x hx l ...dx J 

(18) 

The modified Bessel function Kj-i (A n) can be expressed as multiple deriva- 
tives of K (Xn) jH2]: 

A V An)(^ = (-l)¥(jL)* Ko( A„). (19) 
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Substituting ((T§|) and 1)18(1 into (fT5)) yields Aj(q,X) as an infinite sum over 
the modified Bessel function Ko(Xn) : 



A,(g,A) = 2^rf(-l)^-(-— ) 2 £ K (A n) + ^ ( g , A) . (20) 



dXKXdX 



71=1 



The infinite sum over the modified Bessel function 1£q(A n.) has the following 
series expansion |52| : 

^^o(An) = -{C + ln(A/4vr)} + ^- + vr^ ^=__-- . 

~ 2 2A 1 \ V A 2 + 4 m 2 7r 2 2 run J 

n=l m=l v v ' 

(21) 

By substituting l[2T|) into ((20]) we obtain Aj(g, A) as an analytic expression 
plus the remainder Rj(q, A): 

AJq, A) = — 1 2> vr^ r(2±i) + -4t 2 j+1 VT2 r(2#) 
3W ' ; A^ v 2 y AJ+ 1 v 2 y (22) 

+ A2J+ 2 r(^)7r^-(A) + fl,-(9,A) 
00 j 

where Xj W = Yl ~ 9 , 9 9 ,i±l ■ ( 23 ) 

m=l (A +4™^ 7T^) 2 

To obtain the regularized vacuum energy Ep S9 (q, A) given by (jHJ we need to 
evaluate the derivative of A: 

9aA;(<Z, A) = 2* ^ r(i±i) - ^±1 2 i+ 1 vri r(2± 2 ) 

+ 2^ 2 r(^),r* Xj(A) + A r( i+2 )7r i 5aX . (a) (24) 
+ d x Rj(q,X). 

We now take the limit as A — > in (]24[). Note that the first two terms in 
(j24j) are divergent in this limit and represent the infinite continuum energy of 
surface and volume terms respectively. The Casimir energy is the difference 
between the discrete and continuum case and therefore these two terms need 
to be subtracted out. We therefore define 

d x A f / mte (q,X) = r(i±2)7ri Xi(A)+A 2^' +2 T{^)^ d\ Xj (A) +d\Rj (q,X). 

(25) 

The above terms in the limit A = are 

lim Y i (A) = (2vr)-^ 2 C(j + 2) ; lim «9 A X , (A) = (26) 

A — A — *U 
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and we define Rj(q) as 



Rj(q) = ! im dxRj^X). (27) 

Substituting (|26j) and ()27|) into ()25|) we obtain the compact form 

lim d\Aj inite (q, A) = r(2±2) vr^ CO" + 2) + ■ (28) 

Rj(q) is the multidimensional remainder which contributes to the Casimir 
energy. This is evaluated in Appendix A and the result is: 



£3 /j-„\ A A n— Kj+iQirny/ei + ••• + P m 



w^E^'C;' E E 

m=l v 7 n=l4 m =i 



(29) 

Note that i?j ((/) is zero for j = q. The above expression ()29|) for the remain- 
der is highly convenient. First, it converges rapidly. The Bessel functions 
decrease rapidly and therefore only the very first few numbers in each sum 
are needed to reach high accuracy. Secondly, clever algorithms for Bessel 
functions are well incorporated in many software packages making numerical 
computation of the remainder easy and accurate. The finite part of © in 
the limit A = yields the Casimir energy for the periodic case: 

EpM = flLi E JmfcAf Ue (q,X) 

= -^i n i« e r (^) ^ cu + 2) + rm 

i=l j=q 

with Rj(q) given by (|29j) . Equation ()3U|) is the Casimir energy of a massless 
scalar field moving with velocity v in a d-dimensional rectangular box with 
periodic boundary conditions where d — q sides have length L and q sides 
have much larger lengths. Note the convenient break-up into two terms: 
a finite analytical formula over the well-known Riemann zeta and gamma 
functions plus a remainder. Since Rj(q) is zero for j = q, the sum for the 
remainder starts at j = q + 1 and is therefore non-zero only if d > q + 2 i.e. 
non-zero only if there is at least two small dimensions on top of the q large 
dimensions. 

We can now readily express the Casimir energies for the Neumann and 
Dirichlet cases as sums over the periodic ones. In ©, the sums for the 
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periodic case start at — oo, while for Neumann and Dirichlet cases they start 
at and 1 respectively. We can express the sums from or 1 to oo in terms of 
sums from — oo to oo. The functions being summed have the property f(n) = 
/(— n). We therefore have the relation f(n) = ^ J2°°oo f( n ) + \ /(^) 
which can be expressed as an operator — > \\Y^ o +-*-)■ Applying the 
operator d — q times yields: 



i=l JU i=l J ~°° -oo 

q d—q , „ oo 

i=l m=l — oo 

= 2-^ d f / ( d - q )E p ( q , q + m) 



m=l 



(31) 



Substituting (|30j) into (|31|) yields the Neumann Casimir energy: 



i=i 

-1 d-q 



E E (^(rm-r^CW + !■)+*,(«) 



(32) 



m 

J=9 m=j-q+l 



For the Dirichlet case, Y^^° /(n) = ^ X^oo /( n ) ~~ 2 /(^) ano - we obtain 

d-q 

d — q 



E D { q ,d)=2- d - 1 Y J (-^) d ' q+m { m q )E p (q,q + m). (33) 

m=l 

Substituting (|3*U)) into (|3*S|) yields the Dirichlet Casimir energy: 

E D (q,d)=2- d - 1 ^ T f[L i 

i=i 

EC-i^CT!: 1 ) (r(^)vr^ CO' + 2) + 



(34) 



A special case is that of Dirichlet conditions for parallel plates where all sides 
except one are large i.e. q = d— 1. Rj(q) is then zero and only j = d — 1 is 
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summed: 



E\\ (d) = -2-"- 1 ^ J] L * r (^) C(d + 1) (35) 

The Casimir pressure for the parallel plates is then: 

where V is the volume L Y[i=i °f the parallel plates and we have re- 
inserted ?i. The result ()36j) is in agreement with the higher-dimensional 
parallel plate cut-off calculation of 36] if we set v and L to unity. For three 
dimensions we set d = 3 and obtain: 

fl|W = -s£ (37) 

where we used the fact that £(4) = 7r 4 /90. This result is in agreement with 
the Casimir calculation for quantum fluctuations in a dilute Bose-Einstein 
condensate at zero temperature that was recently carried out by [T51 ITU]. 
As previously mentioned, though the BEC has a non-linear dispersion re- 
lation the Casimir energy only picks out the low frequency part since the 
higher frequencies act as a continuum. The low frequency part is linear and 
the dispersion relation is equivalent to that of a massless Klein-Gordon field 
with speed of light replaced by speed of sound. The pressure in (|37|) is neg- 
ative implying attraction and decreases to the fourth power of the distance 
as in the electromagnetic case. In fact, the classic electromagnetic result 
— 5io TF ^ or P arane l-pl a tes can be obtained by multiplying (|3*7|) by 2 for two 
polarizations and setting v equal to 1 for the speed of light. 

Equations (|3*0|) . (|3*2*)l and (|3*1)) for the Casimir energies contain products of 
the large dimensions Lj which can be arbitrarily large. It is of more physical 
interest to obtain the energy densities e which depend on L only. Dividing 
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the Casimir energies by the volume V = L d q Y\i=i Li yields 



d-l 



L d+i 



TT V 



-j-4 
TT 2 



C(j + 2)+R j (q) 



3=1 



(2L) d + l 



TT V 



3=1 m=j-q+l 




(2L) <*+l 



TT V 




3=1 



(38) 



The three equations in (|38|) are our final results for the periodic, Neumann 
and Dirichlet Casimir energy densities for massless scalar fields moving with 
wave velocity v in a d-dimensional rectangular cavity where d — q sides 
have equal length L and q sides have much larger length. The expressions 
contain a dominant finite analytical part plus a fast-converging remainder 
Rj{q) given by (|29j) . General formulas for arbitrary lengths are obtained in 
appendix B. 

3 Physical interpretation of Casimir energy for- 
mulas 

The Casimir energy formula (|30j) for periodic boundary conditions and (|32jl 
and (|34j) for Neumann and Dirichlet conditions respectively have a clear 
physical picture or interpretation. Excluding the remainder, the formulas 
can be viewed as the energy needed to set up the parallel plates from which 
the rectangular cavity is constructed. For example, consider the case d = 3 
and q = corresponding to a cube (hypertorus for periodic) with sides of 
length L. The cube is built out of three sets of parallel plates. In (|3Ujl this 
corresponds to summing the term r(^p) tt^~ £(j + 2) for j = 0, 1 and 2. 
To build the cube, one begins by placing two plates a distance L apart. This 
corresponds to j = 2. Adding two more plates corresponds to j = 1 and the 
last two plates completes the cube and corresponds to j = 0. We now show 
mathematically that the Casimir energy is the sum of parallel plate energies 
plus a remainder. Consider periodic boundary conditions. The energy for 
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parallel plates defined by letting q = d — 1 in (j3"U|) is: 

E Pn( d ) = -l^fl L i n^)^C(d+l). (39) 

i=l 

i?j (g) is zero for parallel plates and this is why it is not present in (|39j) . The 
parallel plate energy in j + 1 dimensions is 

^I|0-+1) = -T^]I^ r(^)vr^C(j + 2). (40) 

8=1 

In ()3fljl . j > g. Therefore the first g products in nf =1 are large and the 
rest are equal to L so that the above product Yll =i Li can be replaced by 
L J - q ntl Li yielding 

^n0' + i) = -^f[^ r(i±2 )7r ^ c(j + 2) . (41) 

1=1 

Substituting JUJ) in (JHDJ) yields: 

d-l 

E p (q, d)=J2( E P]] U + 1) + ) ( 42 ) 

As can be seen, the Casimir energy in a d-dimensional space with q large 
dimensions is the sum of parallel plates immersed in different dimensions 
plus a remainder. When building the rectangular cavity out of successive 
parallel plates, the first parallel plates have d — 1 large dimensions, the 
second have d — 2 large dimensions and so on until the last set which has 
q large dimensions. In short, the d — q dimensional resonator is the sum of 
one-dimensional resonators each immersed in a different dimension ranging 
from q + 1 to d — 1. 

What is the physical interpretation for the remainder? The energy for par- 
allel plates are by definition those for isolated plates in vacuum. However, 
to construct the rectangular cavity, one adds plates to other plates already 
present. To clarify this difference consider two scenarios. Scenario I: plates 
are brought together in vacuum in a two dimensional space. This leaves one 
dimension which is large. Scenario II: consider a three dimensional space 
where there is already a pair of parallel plates. Now add another pair of 
plates. This leaves one dimension which is large as in scenario I. The main 
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point is this: the energy in scenario II for adding the second set of plates is 
almost but not exactly equal to the energy of the plates in scenario I. The 
reason is that in scenario II there is also an interaction energy due to the 
presence of the other plates. The remainder term is therefore an 'interac- 
tion' or potential energy arising from the nonlinearity of the energy when 
waves moving along different directions are added. By interaction energy we 
do not mean that there is a Feynman diagram where scalar fields meet at a 
vertex. That would be a nonlinear theory like X(p . What we have here is a 
linear theory and the waves obey the superposition principle. However, the 
energy is clearly not linear. This is reminiscent of what occurs in classical 
electrodynamics. In vacuum, the theory is linear and one can add two elec- 
tric field vectors but the energy itself is not linear since it is proportional to 
the square of the electric field. What we usually call the potential energy 
between two static charges q\ and q2 is nothing but the interaction energy 
between the electric field Ei produced by the first charge and the electric 
field E2 produced by the second charge. The energy density is proportional 
to (Ei + E 2 ) 2 = El + E$ + 2 Ei • E 2 and the integration of the cross-term 
2 Ei • E 2 over all space yields the well-known potential energy proportional 
to qiq^/r where r is the distance between the charges. The remainder term 
is similarly a potential energy arising from the nonlinearity of the energy. 

We can now make predictions about the behaviour of the remainder for 
periodic, Neumann and Dirichlet boundary conditions. We predict the fol- 
lowing: 

• percentage wise, the periodic case will have the highest remainder, the 
Dirichlet case the smallest, and Neumann in between 

• the remainder grows with the space dimension for the periodic and 
Neumann cases but actually decreases for the Dirichlet case 

Let us see how we can make such predictions. The Casimir energy is the dif- 
ference between discrete and continuum modes. As the frequency increases 
the discrete approaches the continuum. Therefore the Casimir energy picks 
out the low frequency or low energy behaviour. Moreover, the lower the en- 
ergy, the more nonlinear is the change in energy. Higher energies are closer 
to the continuum and changes are more linear. The minimum energy mode 
for the periodic and Neumann cases is zero (the case when all rtj's are zero). 
For Dirichlet the minimum energy mode occurs when all n^s are equal to 1. 
For concreteness let the space dimension be 5. For periodic and Neumann 
the smallest nonzero energy state occurs when one rii is 1 so that one of 
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five slots is filled with 1 e.g. (0,1,0,0,0) while for Dirichlet the minimum 
energy starts at (1,1,1,1,1). Now add 1 to both cases (creating states with 
two l's like (0,1,0,0,1) and states like (1,2,1,1,1)). The percentage change in 
the energy in the Dirichlet case will not be large because the energy started 
off large. The energy changes almost linearly leading to a small remainder. 
As the dimension increases, the energy for the Dirichlet case starts off even 
higher and the change is even less. For Dirichlet, we therefore predict the 
remainder to be a very small percentage of the energy and that it decreases 
as the space dimension grows. In the periodic and Neumann case, the energy 
starts off low, so the change is a larger percentage of the initial energy and 
therefore more nonlinear than in the Dirichlet case. This effect is greatly 
accentuated by the fact that are many more low-energy combinations for the 
Neumann and periodic case compared to the Dirichlet case. For example, 
there are 5 ways to place 2 in (1,2,1,1,1) but there are 10 ways to arrange the 
two l's in (0,1,0,0,1). The remainder will therefore be considerably larger in 
the Neumann and periodic case. Moreover, the remainder for periodic and 
Neumann cases will grow as the dimension increases because as the number 
of zeros increases there are simply more possible low-energy combinations 
and this increases the nonlinear effect. Finally, the periodic case has the 
largest remainder of all the cases because negative n's are allowed, so that 
in our state (0,1,0,0,1) one can also have combinations with —1 leading to 
considerably more low-energy contributions than in the Neumann case. Our 
numerical results confirm all these trends. 



4 Epstein zeta in odd dimensions as products of 
one-dimensional sums plus remainder 

When applied to a rectangular geometry, the zeta function regularization 
technique via analytical continuation yields quickly a finite expression for 
the Casimir energy in terms of homogeneous Epstein zeta functions. The 
subtraction of two infinities does not explicitly appear anywhere in the pro- 
cess. This is a great advantage over the cut-off technique. We use zeta 
function regularization here to obtain quickly an expression for the Casimir 
energy in terms of Epstein zeta functions for the periodic case. Our main 
goal however is to express the homogeneous Epstein zeta function for 3,5 
and 7 dimensions in terms of products of one-dimensional sums plus a small 
remainder. Readers interested in getting a deeper understanding of the 
zeta regularization technique as well as other techniques such as heat-kernel 
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methods are referred to the following books |551 1561 I57j . A sample of 
older and more recent articles where these techniques are applied in var- 
ious contexts ranging from gravitation to condensed-matter can be found 
in [HHl EHl EOl EH EH E31 EH ESI EH1 EH EHl IMl HOI HU For the case of 
rectangular cavities in arbitrary d dimensions these techniques have been 
applied in [HJ |I3J EU . 

Though one can compute a finite numerical result, extra work must be done 
to express the Epstein zeta function in a compact analytical form. Define 
the Epstein zeta function Za(ai, . . . , a d ; s) as: 

oo 

Z d (a 1 , . . . ,a d ;s) = ^ [{aiuif ^ h (a d n d ) 2 ] 3 (43) 

Tli = — OO 

i=l,...,d 

where the prime excludes the case where all n's are zero and absolute con- 
vergence requires Re s > d/2. Our definition differs from the standard one 
by a factor of 2 in the "power i.e. we have —s instead of —s/2. This defi- 
nition is chosen as it simplifies our final expressions. We focus on the case 
of the hypercube, where all the a's are equal and can be pulled out of the 
sum in (|43[) (for simplicity we set them to unity). This yields the homoge- 
neous Epstein zeta function Z d (s). The vacuum energy in d dimensions for 
periodic boundary conditions is trivial to write in terms of Z d {s): 

oo 

E; ac (o,d) = ^ £ („? + ... + „2)V a 

£w ( 44 ) 

= ^Z d (-l/2). 

Now Z d (—l/2) is formally infinite if (|43|) is applied in a straightforward 
fashion. It therefore requires regularization. The keystone of the zeta regu- 
larization technique is analytic continuation and the existence of a reflection 
formula. Like the Riemann zeta function, the Epstein zeta function has an 
integral representation which yields an analytic continuation over the entire 
complex plane except for a pole at s = d/2. The representation leads to the 
following functional relation or reflection formula: 

tt- s T(s) Z d (s) = ir s - d ' 2 T(d/2 - s) Z d (d/2 - s) . (45) 

We therefore obtain that 

Z d (-l/2) = -0.5 Z d (^) r(^±i) vr^ (46) 
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and reduces to the Casimir energy 



E p (0,d) 



— 7TV 



(47) 



2L 



Clearly, the Casimir energy is finite since Z d (— ^— ) converges. The reader 
should appreciate just how quickly the zeta function technique yields this 
result. 

The homogeneous Epstein zeta function Z d {s) can be expressed in terms of 
sums over the arithmetical function r^{n) which is the number of represen- 
tations of an integer n as a sum of d squares without regard to sign or order: 



We can therefore use results from number theory on rd(n) to obtain directly 
formulas for the Epstein zeta function. It turns out that formulas for rd(n) 
which are not complicated exist in 2, 4, 6 and 8 dimensions and these can be 
used to obtain the Epstein zeta function (|48|) as products of one dimensional 
sums with no remainder. The formula for dimension 1 is trivial (by definition 
a Riemann zeta function) but formulas for 3, 5 and 7 dimensions eluded 
number theorists until a major breakthrough in 2002 when Goro Shimura 
developed a systematic way of finding formulas for r^in) for values of d 
up to 8 jSJ]. Unfortunately, the odd-dimensional formulas are much more 
complicated than the even ones. However, one can develop a technique 
where one obtains excellent analytical expressions plus a small remainder 
for £3(3), Z5 (s) and Zj(s). This technique makes use of number theory 
results in 2,4,6 and 8 dimensions and the Euler-Maclaurin formula to fill 
in the odd-dimensional gaps. The remainder which is explicitly evaluated 
turns out small because the odd cases are derived to a large part from the 
even cases. The most important case is of course Z^(—^) since it relates to 
the realistic three-dimensional Casimir energy. We obtain a nice compact 
analytical expression for Z%(s). The analytical part is so accurate that it 
yields the correct Casimir energy to within a remarkable 0.04% as compared 
to 1.6% from our cut-off formulas. 

We start by stating the number-theoretic formulas for r2(n),r±(n) : r§{n) and 
r$ (n) and the known exact expressions for Z\ , Z2, Z4, Z§ and Z% obtained 



00 



Z d (s)= Y! H + --- + n 2 d ] 



m=— 00 
i=l,.,.,d 



(48) 




n=l 
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from them via Q48JI . We illustrate how to obtain Zd(s) via the number- 
theoretic formulas for r^(n), something that may not be too familiar to 
many Physicists. We choose d = 6 as the example to illustrate as it fills 
a gap in the table quoted in [32j which contains Z\,Z%, Z^ and Z$ but not 
Zq. We then develop the mathematical technique by which we obtain the 
odd-dimensional homogeneous Epstein zeta functions. 



4.1 Exact expressions for even-dimensional Epstein zeta func- 
tion via r^in) 

As mentioned already, the arithmetical function r^(n) is the number of rep- 
resentations of an integer n as the sum of d squares without regard to order 
or sign. The formulas for r^(n) for d = 2, 4, 6 and 8 are known and given by 
(a good history with references can be found in |72| ) : 

r 2 (n) = A^xid) r 4 (ra) = 

d\n d\n 

r 6 (n) = 16 ^ x{d')d 2 -4^ X (d) d 2 r 8 (n) = 16 ^(-l) n+d (49) 

d\n d\n d\n 

where d! = n/d and x(d) is the primitive Dirichlet character modulo 4 given 

by x{d) = if d is even and x(d) = ( — l) -5 " if ^ is odd. We now evaluate 
Zd(s) for d = 6: 



Ms) = E 



r 6 (n) 



n=l 



<2'=oddd=l v ; d=oddp=l v ^' 

OO OO / - l \ rY1 oo oo 



16VV i-ir (-ir ( 5 °) 

d tt ( 2 - + ^ 2 P ^ (2- + I)- 2 P s 

^ (2m + l) s ^ d s ~ 2 ^ (2m + l) s ~ 2 ^ p s 

m=0 y ' d=l m=0 y > p=l r 

16/3(sK(s-2)-4/?( S -2)C(s) 



where f3(s) and £(s) are the Dirichlet beta and Riemann zeta function re- 
spectively defined by /3(a) = E^°= (- 1 ) n /( 2n + X ) S and C(«) = E^°=i V" 8 - 
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We have illustrated how knowledge of the arithmetical function r^n) leads 
to an exact and simple representation for the Epstein zeta function Z§{n) 
as a product of the one-dimensional sums (5{s) and C( s )- The other Epstein 
zeta functions can be obtained in a similar fashion. We state them below 
together with Z 6 (s) 

Z 1 (s) = 2((2s) 
Z 2 ( S ) = AC(s)(3(s) 

Z 4 ( S ) = 8C(s)(( S -l)(l-A 1 - s ) (51) 
Z 6 (s) = 16/3(8) C(s-2) - 4/3( S -2) C(s) 
Z 8 ( S ) = 16C( S )C(s-3)(l-2 1 - s +4 2 - s ). 



4.2 Analytical expressions for Epstein zeta function Z^s) in 
3,5 and 7 dimensions 

As already mentioned, the formulas for rd(n) for d = 3,5 and 7 are much 
more complicated than the even ones and it is not easy to use them to obtain 
analytical formulas for Z^^Z*, and Zy. We therefore develop a separate 
technique to find such expressions. The Epstein zeta function Zd(s) defined 
in (|48|) contains d sums which begin at —00. It is convenient to define 
another function Pk(s) as k sums which start at 1: 

00 

Pk(s)= H + --- + nl]- s . (52) 

m=i 

i=l,...,k 

We can express Pfc(s) as sums over Z m (s): 

k , 

pus) = T,(-v m+k2 ~ k {!j ( 53 ) 

m=l 

Similarly, we can express Zd(s) as sums over Pk(s): 

d , 

Zd(s)=^2( k )2 k P k (s). (54) 
fc=i 

It is instructive to map out the main idea or process behind the technique 
we will use. Consider the example of wanting to find expressions for Z3. 
From (|54[). you would need to know P%, P2 and P3. You can find P\ and P2 
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from (|53|) since analytical expressions for Z\ and Z2 are known. However, 
you do not know P3. At this point, you use the Euler-Maclaurin formula to 
express P3 in terms of P2 plus a remainder. Again, you know P2 in terms 
of Z\ and Z2, so that you can finally express Z% in terms of Z±, Z2 and a 
remainder and hence as an analytical part plus a remainder. The process 
can be continued to find expressions for Z5 and Zj (and even Zg if one 
wants to but the expression becomes cumbersome). We now develop the 
mathematical technique and obtain our main equation. Z d (s) given by (|54[) 
can be expanded as 

Z d (s) = Y J ( k )2 k Pk(s) + 2 d P d (s) 

fe=1 (55) 
d-i , 



k=l 



where (|53*)) was used. We now express P d (s) in terms of P d -i(s) plus a 
remainder via the Euler-Maclaurin formula 0: 



00 

Pd(s)= i n i + ---+ n d] 

nj=l 
i=l,...,d 

OO /»rvi 7 

ax 



}Z J o ( x 2 +n 2) S 2n 2, ( 56 ) 

1o 1 



i=l,...,d-l 

d 2 



- 00 „i 



2 ^qJo dx 2 ({x + v) 2 + n 2 ) 



where 



n 2 = n\ + • • • + n\_ x . (57) 
The first integral in ()56|) can readily be evaluated: 



dx _ 1 T(s-i) 



7T 



(x 2 + n 2 ) s n 2 *- 1 T(s) 
a(s) 1 



(58) 



n 2s 



where a(s) is defined by 



_ 0Fr( s -|) 

= wpj ' ( 59 ) 
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Inserting Q58JI into yields 

a(s) 



Pd{s) = E 



t=l,...,d-l 



2n 



2s-l 



2n 2s 



~E/ 



(60) 



c) 2 



9x 2 ((x + z^) 2 + n 2 ) 



■ dx . 



00 j 

By definition V —5- = P^_i(s). Therefore 

n*=l ™ 
i=l,...,<Z-l 



ft« = ^ ft-i(« - I) - \ Pd-i(s) + R d (s) 



where Rd(s) is the remainder defined by 

Rd(s) = 



00 1 00 „! o2 



I 



dx 2 ((x + i^) 2 + n 2 ) 



(61) 



■ . (62) 



ni ,...,n ( j_ 1 =l i^=0 

The remainder Rd(s) is worked out in appendix C and the result is 



Us) = E E 



00 00 o / /)\ s— 1/2 

Z / 7T C N 



ni,...,n d _i=l f=l 



'7T \ n 



r(l-s) sin(7rs)ir s ^ /2 (27r^n) (63) 



where n is given by (|57|). We now evaluate the term 2 rf P^(s) occurring in 
(j55|) via (|SH and lfK5|l : 



2 d P d ( S ) = 2^ a(s) P d _i(s - |) - 2 d - 1 P d _!( S ) + 2 d P d (s) 



m=l 



[a(s) Z m {s - \) - Z m {s)] + 2 d R d (s) . 

(64) 

Substituting into we obtain our main equation: 



^) = E(- 1 ) d+m_1 



m=l 



d-l 



in 



a(s) Z m (s - i) + 



d-l 
m — 1 



Z m (s) 



+2 d R d (s) 
(65) 

where Rd(s) is the remainder given by (|d^|) . Equation (|63|) expresses as 
sums over Zj's from 1 to d— 1 plus a remainder. We are now in a position to 
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obtain expressions for Z^,Z^ and Z-j as products of one-dimensional sums 
plus a remainder by using our main equation (|65|) together with the ana- 
lytical expressions for Z\, Z 2 , Z^, Z§ and Z$ given in (|S*Tj) . We begin with 
Z 3 (s). Applying equation (ffiSl yields 

Z 3 (s) = a(s) [-2 Z 1 {s -\) + Z 2 {s- ±)] - Z x {s) + 2 Z 2 (s) + 2 3 R 3 (s) . (66) 

We now substitute the analytical expressions for Z\ and Z 2 given in (|51|) 
and obtain our final expression for Z 3 : 

Z 3 (s) = 4 a(s) C( S -i)/3(s-i) -4 a(a) C(2s-1) +8 CW(s) -2 C(2s) +8 R 3 (s) . 

(67) 

This is a compact analytical result for the important three-dimensional case. 
The only remainder is 8R 3 (s) and the rest includes four analytical terms, 
each expressed in terms of simple one-dimensional sums and gamma func- 
tions. Later we will see that the analytical part yields numerically the correct 
Casimir energy to within 0.04%! We now evaluate Z^(s). Using again the 
main equation (|65j) we obtain: 

Z 5 (s) = a(s) l-AZ^s-l) + 6Z 2 (s - ±) - AZ 3 (s - ±) + Z 4 ( S - |)] 

c 68 
-Z 1 (s) + AZ 2 {s)-6Z 3 {s) + 4Z i (s) + 2 5 R 5 (s). 

Substituting the analytical expressions for Z\,Z 2 and Z 4 given in (j51j) and 
Z 3 from (|67|) into (|68|) one obtains the final expression for Z^{s): 

Z 5 (s) = 10 C(2s) - 32 ((s)p(s) + 32 ((s) C(s - 1) (1 - 4 1 -*) 

+ 8a(s) [3C(2s - 1) - 4C(s-±)/3( S -±) + COO C(*-f)(l " 2 3 " 2s )] 
- 16a(s)a(s-|)(C(s- l)/3(s - 1) - C(2s - 2)) + RZ 5 (s) 

(69) 

where the remainder RZ 5 (s) = -32 a(s) tf 3 (s-±)-48 R 3 (s) + 32 i? 5 (s). The 
expression for Z-j(s) is: 

Z 7 (s) = a(s) [-6Z x (s-±) + 15 Z 2 {s - \) - 20Z 3 (s - ±) 

+15Z 4 (s-±)-6Z 5 (s-i) + Z 6 (s-±)] -Zi(a) + 6Z 2 (s) (70) 
- 15 Z 3 (s) + 20Z 4 (s) - 15Z 5 (s) + 6Z 6 (s) + 2 7 i? 7 (s) 

where Z\ , Z 2 , Z 4 and Z 6 are given by (|5*Tj) , Z 3 by (|6*7)) and Z 5 by l|E9*|) . It 
would be cumbersome to write out the analytical terms for Zj as we did 
for Z 3 and Z5. For calculations, one simply evaluates the necessary Z's 
and substitutes them in (|7()j) . This ends our results for the odd-dimensional 
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Epstein zeta functions. One could have continued and obtained expressions 
for Zg(s) but this is no longer interesting as the expressions become way 
too long. We now state and discuss the numerical results for the Casimir 
energy. 

5 Numerical results and discussion 

Table 1 contains the numerical results for the Casimir energy density for 
periodic (e p ), Dirichlet(eD) and Neumann(eTv) for q large dimensions and 
d — q dimensions of equal length L. This is calculated using the formulas in 
(|38|) and the equation (|29[) for the remainder Rj(q) (v and L are assumed 
to be unity). We state the analytical and remainder contribution separately 
and calculate their sum to obtain the Casimir energy density. For dimensions 
up to d = 5, we include all values of q. For higher dimensions up to d = 10 
we only state q = 0. For numerical results for the case where one has 
arbitrary lengths the reader is referred to |431 1441 145j . The formulas derived 
in appendix B are actually very well suited for such a numerical study but 
length limitations restrict us here. 

Table 1 shows that the absolute value of the Casimir energy density for the 
periodic case is the largest, followed by the Neumann and Dirichlet. Note 
that the sign in the Dirichlet case alternates in two fashions: for a given q, 
it alternates as the dimension d changes and it also alternates as q changes 
for a given d. The Casimir energy densities agree with a few exceptions 
with results obtained by computing the Epstein zeta function and quoted 
in the table in |l2]. For periodic boundary conditions, results for d = p 
(corresponding to q = in our case) are close to our values but do not fully 
agree. For d = 2 the values agree but for d = 3 they obtain —0.81 while 
we obtain —0.838. For d = 4, they obtain —0.85 while we obtain —0.932 
and for d = 5 they obtain —0.95 while we obtain —1.022. Which values 
are correct? Table 3 contains an independent determination of the Casimir 
energy density for the case q = for periodic boundary conditions. The 
values in Table 3 for d = 3, 4 and 5 are -0.837537, -0.932077 and -1.02283 
respectively and these values are in agreement with our results. Therefore, 
in the few places where our results differ from [121 > our numerical values can 
be considered correct. Some numerical results are also quoted for Dirichlet 
boundary conditions in |431 144 j where Epstein zeta functions were also used. 
In the column u = corresponds to our q = and are in agreement. In 
|44j where D is the spacetime dimension i.e. D = d + 1, their first column 
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corresponds to our d — q = 2 results and are in agreement. 

In Table 2 the percentage of the Casimir energy which is a remainder is 
quoted for the different boundary conditions as a function of the dimension 
d (for simplicity, we quote the hypercube case q = but the same trend is 
followed by all q values). Table 2 confirms the predictions made in section 3. 
Moving down the table, as the dimension increases, the percentage decreases 
for Dirichlet but increases for Neumann and periodic as predicted in section 
3. Moving horizontally across the table the percentage is lowest for Dirichlet 
and largest for periodic with Neumann in between, again as predicted in 
section 3 (with the only exception being d = 2 due to the limited low- 
energy permutations in the periodic and Neumann case and the fact that 
the Dirichlet starts off at a low energy unlike higher dimensions). 

Note how small is the percentage remainder. Only at the highest dimensions 
is the percentage high and this mostly for the periodic case. The percentage 
remainder is negligible for the Dirichlet case and the analytical formulas are 
all we need. The Neumann case has a very low remainder at low dimensions. 
At d = 4 it has less than a 1% remainder so that the analytical formulas 
are simply excellent at lower dimensions. Even the periodic case at d = 3 
has only a 1.6% remainder but the remainder grows rapidly with dimension 
compared to the other two cases. 

Table 3 contains the Casimir energy for the periodic case at q = for values 
of d ranging from 2 to 8 calculated via the expressions for the homogeneous 
Epstein zeta function Z^s) (again v and L are assumed to be unity). Our 
aim here was not to make a complete table of Casimir values using the Ep- 
stein zeta function. This has already been successfully done in [12] • The goal 
was mainly to calculate the analytical and remainder terms for the homoge- 
neous Epstein zeta function in 3, 5 and 7 dimensions. For even dimensions, 
the expressions are calculated via (|51|) where there is no remainder. For 
the odd cases of 3, 5 and 7 dimensions they are calculated via our derived 
expressions (IBT)) . (|6U|) and ((70)) and (|6^)l for the remainder Rd(s). Note how 
close are the derived Epstein zeta analytical results to the actual Casimir 
energy and hence the small remainder percentage wise. The analytical ex- 
pressions (|6TJ) . (|69j) and (fTflj) we derived for the Epstein zeta are limited to 
a few dimensions but are exceptionally accurate. As already stated, for the 
realistic three-dimensional case, the remainder is only a remarkable 0.04% 
of the Casimir energy. As one can see, the remainder for these few cases is 
smaller than the remainder from our cut-off technique. The reason is due 
to the fact that the odd-dimensional cases are derived from the even ones 
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Table 1: analytical and remainder contributions to Casimir energy density 
for periodic, Dirichlet and Neumann boundary conditions 



(q,d) 


Periodic 


Dirichlet 


Neumann 




Ep 


analytical 


remainder 


Ed 


analytical 


remainder 


En 


analytical 


remainder 


0,2 
1,2 


-0.718873 
-0.191313 


-0.7149121 
-0.1913133 


-0.0039607 



0.041041 
-0.023914 


0.0415357 
-0.0239142 


-0.0004951 



-0.220759 
-0.023914 


-0.2202637 
-0.0239142 


-0.0004951 



0,3 
1,3 
2,3 


-0.837537 
-0.305322 
-0.109662 


-0.8245743 
-0.3009756 
-0.1096623 


-0.0129625 
-0.0043463 



-0.015732 
0.004832 
-0.006854 


-0.0156650 
0.0051032 
-0.0068539 


-0.0000675 
-0.0027164 



-0.285309 
-0.042997 
-0.006854 


-0.2837567 
-0.0427251 
-0.0068539 


-0.0015528 
-0.0002716 



0,4 
1 4 

2,4 
3,4 


-0.932077 
-0.394299 
-0.193407 
-0.078797 


(\ nmn i a 

-\).y\)jj 1 14 

-0.3797726 
-0.1884593 
-0.0787971 


-U.UZo 

-0.0145260 
-0.0049473 



0.006226 
-0.001634 

0.000810 
-0.002462 


U.UU0Z4J J 

-0.0015871 
0.0009645 
-0.0024624 


-0.0000465 
-0.0001546 



-0.334058 
-0.058881 
-0.012898 
-0.002462 


-Kj.jjxj 1 VoU 

-0.0580200 
-0.0127432 
-0.0024624 


-U.UUjZuUU 

-0.0008614 
-0.0001546 



0,5 
1,5 
2,5 
3,5 
4,5 


-1.02283 

-0.478283 

-0.270975 

-0.150257 

-0.065622 


-0.9689932 
-0.4453944 
-0.2540811 
-0.1444188 
-0.0656218 


-0.0538395 
-0.0328889 
-0.0168939 
-0.0058383 



-0.002611 
0.000504 
-0.000308 
0.000115 
-0.001025 


-0.0026055 
0.0005171 

-0.0002764 
0.0002059 

-0.0010253 


-0.0000052 
-0.0000135 
-0.0000321 
-0.0000912 



-0.372895 
-0.072698 
-0.018440 
-.0004810 
-0.001025 


-0.3671673 
-0.0708686 
-0.0017945 
-0.0047190 
-0.0010253 


-0.0057281 
-0.0018292 
-0.0004959 
-0.0000912 



0,6 


-1.12249 


-1.029970 


-0.092517 


0.001114 


0.0011158 


-0.0000017 


-0.405594 


-0.3964942 


-0.0091000 


0,7 


-1.24313 


-1.091817 


-0.151318 


-0.000489 


-0.0004884 


-0.0000006 


-0.434680 


-0.4211207 


-0.0135591 


0,8 


-1.40015 


-1.159323 


-0.240830 


0.000217 


0.0002170 


-0.0000002 


-0.461950 


-0.4426076 


-0.0193419 


0,9 


-1.61621 


-1.237827 


-0.378385 


-0.000098 


-0.0000977 


-0.0000001 


-0.488792 


-0.4620371 


-0.0267547 


0,10 


-1.92725 


-1.334376 


-0.5928761 


0.000044 


0.0000444 


0.0000000 


-0.516394 


-0.4801973 


-0.03619642 



which contain no remainder. 
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Table 2: percentage of Casimir energy which is remainder (case q=0) 



d 


Dirichlet 


Neumann 


Periodic 


2 


1.2 


0.22 


0.55 


3 


0.42 


0.54 


1.6 


4 


0.31 


0.98 


3.1 


5 


0.20 


1.5 


5.3 


6 


0.15 


2.2 


8.2 


7 


0.12 


3.1 


12 


8 


0.09 


4.2 


17 


9 


0.10 


5.5 


23 


10 


0.00 


7.0 
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Table 3: Epstein-zeta function and comparison of remainder with cut-off 



d 

(q=0) 


Casimir 
energy density 

(periodic) 


Epstein Zeta 
Analytical Remainder 


% remainder 
Epstein zeta Cut-Off 


2 


-0.718873 


-0.718873 








0.6 


3 


-0.837537 


-0.8372276 


-0.0003928 


0.04 


1.6 


4 


-0.932077 


-0.932077 








3.1 


5 


-1.02283 


-1.025582 


0.0027514 


0.3 


5.3 


6 


-1.12249 


-1.12249 








8.2 


7 


-1.24313 


-1.197224 


-0.0459060 


3.7 


12.2 


8 


-1.40015 


-1.40015 








17.2 
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A Remainder term Rj(q) 



In this appendix we evaluate the remainder term Rj(q) defined by 

RM = 

where Rj(q,X) is given by (|17l) i.e. 



Rj{q)= lim dxRjfaX) (71) 

A — *0 



j-q oo . _ rOQ „ i rn oo ^ 2 

%W)=EEHr( , „') !, - tI | / IIE B ^.>'" 

m =l ra =l u u i=li/j=0 

e- A ^ a +^+^) a +-+(*™+^) a +^+-+^- m ,,,,...,/,,,.,,„.... dy ._ m . 

(72) 

There are m integrals from to 1 over the x's and j — m integrals from 
to oo over the y's. Our goal is to simplify (|72|) as much as possible and put 
it in a compact form useful for computations. In the end, the result is that 
(|72[) can conveniently be reduced to sums over Bessel functions. The first 
step is to convert the multiple integrals over the y's to a single integral by 
using spherical coordinates: 

j — m 

r 2 = y\ + ---+ y). m ;d yi ... d Vj . m = 2"^' +1 r^ m ~ l dr . (73) 

Rj(q, A) is then reduced to 

J'-g 00 /n-n\ tt^t- r°° r 1 m 00 Pi 2 

e -A Vn 2 +(.i+.i) 2 +-+(x m +%) 2 +r 2 ^ r J-m-l dr _ 

(74) 

We now turn to the x- integrals from to 1. Note that x + v is continuous 
and runs from to oo. It is therefore convenient to drop the sum over v, 
replace x + v by x and integrate from to oo instead of to 1. This is valid 
as long as the Bernoulli function B-2(x) is replaced by B2(x— [x]) where [x] is 
the greatest integer less than or equal to x. This ensures that the Bernoulli 
function is periodic with period 1 while x runs to infinity. Moreover, B-2(0) = 
-62(1) so that Bi{x — [x]) is not only periodic but continuous. A fourier 
expansion of B 2 {x) = x 2 -x + 1/6 can readily be obtained and is given by 

9 -, /„ v^COs(2 7T^x) 

x 2 -x + l/6 = Y J 2 o (75) 



7T^ 
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The right hand side of (|75jl is a continuous periodic function valid for all x. 
It is equal to the left hand side only for < x < 1 but equal to B2 (x — [x] ) 
over the entire region of integration < x < 00. We can therefore make the 
following replacement: 



dxi 

(76) 



£ / Mxi) 9 f ^l +U ^ dxi [ cos(2vr^x, 

where /(xj) is the exponential function in (|74|) with ^ omitted i.e. 

/(xj) = e" A y/ n2 + x i+-+ x i+-+ x2 m+r 2 . (77) 
The function / has the following properties: 

lhn tlfM = 0; lim ( ?Jp± = Q ] lim /(x f )=0. (78) 

After integrating by parts twice and using the above properties of /, (|76jl 
reduces to 

f cos(2ir£iXi) — 1}~ — dxi = — 4-7r 2 l 2 f cos(27r^j xi) f(x{) 
Jo dxi Jq 

and V / B 2 (z) ^ j ^ -> -4 V / cos(2^x)/(x) dx . 

(79) 

Substituting ((79*)) into equation ((71)) yields 

j— <3 00 00 . i~ m poo m 

m=lfl=l<i ... m =l 1 I 2 ^ 1/0 i=l 



e -A V™ 2 +^+-+^+ f - 2 d Xl . . . dx m r^-" 1 - 1 dr . 

(80) 

We can reduce the above expression ()8U|) to sums over the modified Bessel 
function Kj+i by applying sequentially the following set of three integrals 

EH: 

\b 



I. J oo cos( 7 x)e- Av/ ^+ ? dx= K_ 1 (b v /\ 2 +-f 2 ) 

V A 2 + 7 2 
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II. J °°(x 2 + b 2 ) T 2 u K^aVx 2 + b 2 ) cos(cx) 

(|)i/2 a*" bh* v (a 2 + c*)^"-* K ±u _i(bVa 2 + c 2 ). 

x 2fl+1 2T(/i + l) 

III. / °° K v (a v^ 2 + x 2 ) — - x2)y/2 dx = a)i+1 tf„_„_i(a z) . 

Integral I is applied once and converts the exponential and one cosine into 
the modified Bessel function K-\ i.e. 

r cos(2 7T h xi) e" A dxi = - A6 K_i ( b J\ 2 + 4tt 2 l\ ) 

VA 2 + 47r 2 q v 

(81) 

where 6 = \/ x 2 , + • • • + + n 2 + r 2 . We now make repeated application 
of integral II for the remaining x's that appear in the definition of b. The 
subscript of the Bessel function is therefore decreased by 1/2 each time. 
Since there are m — 1 x-integrals to perform, and we start with K-i, this 
yields the Bessel function K - m -i i.e. 

2 

/ — TTcos(27rl^) b K_ 1 (bJ\ 2 + 4ir 2 £ 2 1 )dx 2 ...dx m 

JO V A + 47F ^1 i=2 

= + ^ 4 = K^i(2vr7^T^^ + ^ + ... + C) 

vr2-(^+^ 2 + ...+C) 4 

(82) 

We now apply integral III to perform the integration over r i.e. 



7r 9 m + 1 i~ m i±l 



(83) 
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The integrals over x and r appearing in (|5U)) can now be replaced by (|53*|) 
yielding: 



.7 + 1 



m=l v 7 n=l^i..... m =l 



(84) 

Finally, by taking the derivative of Rj(q, A) with respect to A and taking the 
limit as A — > yields our desired final result for the remainder Rj(q): 

Rj(q) = lim d x Rj{q,X) 
A — *0 



(85) 

Our final expression (|85|) for Rj(q) is excellent for numerical calculations 
because it converges very quickly (exponentially fast). The sums to infinity 
are formalities as one can reach an accuracy of 8 to 10 digits by summing 
fewer than 9 numbers in each sum for j up to 10. 



B Casimir energy in rectangular cavities with ar- 
bitrary lengths 

One can generalize the multidimensional cut-off method used in section 2 
to obtain Casimir energy formulas for arbitrary lengths in a d-dimensional 
rectangular cavity. Our analysis will naturally be brief since it follows closely 
that of section 2 and many results from that section can be applied here. The 
best way to read this appendix is therefore to have section 2 and appendix 
A in hand for immediate reference. 

The quantized frequencies u for periodic (p), Neumann (N) and Dirichlet 
(D) conditions are now given by: 

u; p = 2nv(§ + --- + §)^ 

(86) 

where the lengths range from L\ to Lj- The regularized vacuum energy for 
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periodic boundary conditions is then given by a similar form to (jSJ) i.e. 



\ I d J. M ^ °° \ / ^ 

TtH hyr / -Ai TT 



^(A) = -nvdx E e Ld =-™a A (i+ 5Z e V L i 



n;=— 00 «i=- 00 

i=l,...,d 



00 00 , l n l , A 00 oo \n\ n\ 



+ E' E e V 1 ? ^ + •••+ ^' ^e - A \/^ + ^ 

n2=— oon\=— 00 ntj=— 00 n^=— 00 

i=l,...,d-l 



d-l 

-Try Y,dxAj(X) 

3=0 



(87) 



where 

00 00 . I n 2 n\ rij 

A,(A) S E' E ^ ^- 

ra=— 00 rii=— 00 
i=l,..., j 

As in (1121). we obtain via the Euler-Maclaurin formula that 



00 />oo 

E /("<) = / f{x)dx-R. (89) 



n,=— 00 



i? is given by expression Q79|) obtained in appendix A: 

00 />oo 00 ^ />oo 

= -4V cos(2vr^a;)/(a;)(ix = -2 V / cos(2 n I x) f (x) dx 
1=1 Jo e=-oc Jo 

(90) 

where we used f{x) = f(—x) for the function we are considering. Then (|89|) 
reduces to 



53 f(m) = 2 53 / cos(2^£x)/( 

n,=-oo £=-oo 



x) dx (91) 



where £ = is now included. Therefore the j-dimensional sum appearing 
in (|88j) for Aj(A) can be obtained by repeated application of (|9*Tj) . What 
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appears in the regularized energy l|87|) is the derivative d\ Aj(A): 

OO OO 

d x A j (X) = d x J^ 



oo oo , / n 2 . «i . 



n=— oo rii=— oo 
i=l,...,j 



00 00 /-oo * / " 2 1 x l 1 

d x ) 2 J > / cos(2 7T^xi)...cos(2 7r £j Xj) e V j+i 1 3 dx\ . . . dxj 



n=-oo l i= — 00 
1=1,..., j 



= T^OT ( 2J+1 ^E r e-V" 2 ^ 4 -^ . . . ds,- + d\Rj(\) ) 

(92) 

where the sum over all £'s was divided into two cases leading to the two 
terms in the brackets of 1)92(1 . The first term occurs when all £'s are equal to 
zero. The second term is for all other fs and corresponds to the remainder 



dxRj(X) 



00 00 



2 j+1 d x J2 E' / «*(2t*1 7^x1)... cos(2vr^ ^ a;,-) e A \Z n2+x2 i+"+ x 3 ' dxi . . . dxj 

n=ll l =-od' J 3 

i=l,...,j 

(93) 

where the prime over the multiple sum excludes only the case when all ts 
are equal to zero. The multiple integral over j cosines can be obtained 
directly from in appendix A by the following substitutions: m — ► j, 
ti — ► ^ Li/Lj+i and n 2 + r 2 — > n 2 i.e. 

a A i?,(A) = 



00 00 2A(nL J+1 ) — ^. +1 (|^^i^+iL + ( 4Ll) 2 + ... + ( ^ Lj)2) 

E ^ zh 

n=l/ ?: =-oo _ /(AL, + i) 2 , ln T N9 , , f „ T ^\ 4 

(94) 



The Casimir energy is proportional to the finite part of (|92j) as A — > 0. The 
first term in brackets in (|92j) is identical to the derivative of the first term 
in Aj(q, A) given by (|15|). Therefore the result ((28|) from section 2 is directly 
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applicable i.e. 



lim d x Af ™ te (A) = (r(^) tt^F 4 C(i + 2) + fl,) (95) 

where the remainder term is given by 



Ri = lim dxRi(X) 



oo 



2 (n 



j'+l 



E E t^i±iL^ — m K m (^-j(e 1 L l r + ... + (e J L J r ) 

n=lli=-oo n [(£ 1 L 1 ) 2 + --- + (£ j L j ) 2 }~ 2 



i=l,...,j 



(96) 



Our final Casimir energy expression for periodic boundary conditions is then 
given by 



e Pli Ld {d) = -,vY J ^d x Kf^{x) 

3=0 
d-1 



-E7^r( r (^)^CO- + 2) + ^ 

3=0 V 3+1 ' 



-it Li C(3) LiL 2 ir 2 ttLi 7tLiL 2 

6Lx L| 2tt L| 90 + "' 1 L| 2 L| + 

(97) 

where the remainder Rj is given by (|96|) (note that i?j is zero when j = 0). 
Equation (|97j) is a highly compact way to express the Casimir energy for 
arbitrary lengths. As in section 2 it is split into two terms: an analytical part 
and a remainder. The same physical interpretation follows: the analytical 
part is a sum of parallel plate terms. Equation (|97|) is valid for any lengths 
and we know the result should be invariant under a permutation of the 
lengths. However, the two terms separately are not invariant, only their 
sum. We naturally want to label the lengths such that the remainder term 
lives up to its name. This can be accomplished if the largest length is 
labeled Li, the next largest length L 2 , i.e. L\ > L 2 > L3.... Then the 
Bessel function decreases exponentially fast and the remainder is small. If 
q dimensions are large and d — q dimensions have equal length L, Ea. (|97l) 
for E p (d) and Eq. (|HH|) for the remainder Rj reduce to the results of section 
2 i.e. E p (q,d) given by (|3*0|) and Rj(q,d) given by ((2*9")) respectively. 
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The Neumann (N) and Dirichlet (D) cases can be obtained via simple per- 
mutations of the periodic case. The operator relations for Neumann and 
Dirichlet are ^(E-oo+ 1 ) and YT ~> ifX^oo -1 ) respectively. 

Applying the operator d times while keeping each sum distinct because of 
different lengths and multiplying the final result by ^ yields the Neumann 
and Dirichlet energies 

m=l (ki,...,km) 

where the (+) is for Neumann and the (-) for Dirichlet. The sum is over 
all sets (ki, . . . , k m ) with k\ < k 2 < • • • < k m (the fc's run from 1 to d). 
E p (m) is the periodic energy (|97|) replacing d by m and L\ by L\. , L 2 
by L k2 , etc. 



C Remainder term Rd(s) for Epstein-zeta function 

We derive in this appendix a convenient form for the remainder Rd(s) in 
terms of sums of Bessel and gamma functions. We begin with the expression 
for the remainder Rd(s) given by (|62[) : 

w ^ 2 ^ J Q y 1 dx 2 {(x + v) 2 + n 2 ) s K ' 

ni,...,n ( ;_ 1 =l ^=0 x 

where 

n 2 = n? + • • • + n\_ x . (100) 

We now follow similar procedures as those employed in appendix A for 
Rj(q)- To avoid being repetitive, we skim through details already discussed 
in appendix A. 

The term x + v is continuous and runs from to 00. We drop the sum over 
v, replace x + v by x and integrate from to 00 instead of to 1. We replace 
B 2 {x) = x 2 — x + 1/6 by its fourier expansion (|75|l i.e. 



x 2 -x + l/6 = ^ ^ ; . (101) 
We can therefore make the following replacement in (|99|) : 

E r^(x) Pfa + ^ dx^J-s r cos(2n lx)^±dx (102) 
^ 7 <9x ^ t 2 ir 2 Jo dx 
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where f(x) is the function in l|9*9")l with v omitted i.e. 

The function f(x) has the following properties: 

Jim ^JM. = ; lim = ; lim f(x) = . (104) 

x^O OX x^oo OX x— >oo 

With the above properties of /, (|102|) reduces to the same expression (|T§|) 
obtained in appendix A: 

J2 B 2 {x) g x dx ^ / cos(2vr^x)/(x)dx. (105) 

After substituting ()105|) into (|9"9*)) we obtain Rd(s) in the following form: 



Ms)= J2 Y\ I 2 cos(2ttIx) - dx. (106) 

ni,...,n d _i=l £=1 1/0 ^ ' V 

The integral can be expressed in terms of Bessel functions i.e. 

f°° 12 /7rA s ~ 1 / 2 

/ 2 cos(2 7r£x)— ^ ^— = — = r(l— s) sin(-7r s) K s _i/ 2 (2tt £n) I — 

Jo \ x +n z ) s yJ-K 1 \ n J 

(107) 

Our final expression for Rd{s) is then 



2 /7rA 
rn,...,n d _ 1 =l £=1 v 71 " V n / 



s-1/2 



(108) 



where n = J n\ + • • • + . 
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